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Chern insulators are band insulators exhibiting a nonzero Hall conductance but preserving the 
lattice translational symmetry. We conclusively show that a partially filled Chern insulator at 1/3 
filling exhibits a fractional quantum Hall effect and rule out charge-density wave states that have not 
been ruled out by previous studies. By diagonalizing the Hubbard interaction in the flat-band limit 
of these insulators, we show the following: The system is incompressible and has a 3-foId degenerate 
ground state whose momenta can be computed by postulating an generalized Pauli principle with 
no more than 1 particle in 3 consecutive orbitals. The ground state density is constant, and equal 
to 1/3 in momentum space. Excitations of the system are fractional statistics particles whose total 
counting matches that of quasiholes in the Laughlin state based on the same generalized Pauli 
principle. The entanglement spectrum of the state has a clear entanglement gap which seems to 
remain finite in the thermodynamic limit. The levels below the gap exhibit counting identical to 
that of Laughlin 1/3 quasiholes. Both the 3 ground states and excited states exhibit spectral flow 
upon flux insertion. All the properties above disappear in the trivial state of the insulator - both 
the many-body energy gap and the entanglement gap close at the phase transition when the single- 
particle Hamiltonian goes from topologically nontrivial to topologically trivial. These facts clearly 
show that fractional many-body states are possible in topological insulators. 



INTRODUCTION 

The interest in the field of topological insulators has 
exploded in recent years, fueled by the theoretical predic- 
tion fT\ and experimental observation pj H] of topolog- 
ical insulators that preserve the time reversal symmetry. 
These insulators have an insulating bulk but exhibit per- 
fectly metallic edge or surface states vifhose existence is 
required by the nontrivial topology of the bulk electronic 
structure. However, most of the theoretical work on topo- 
logical insulators is essentially of single-particle nature 
and neglects interactions, except for perhaps at random 
phase approximation level. The few interacting states [S] 
that have been proposed to give rise to topological insula- 
tors can be understood from a mean-field perspective 
- their quasiparticle excitations are still electronlike and 
not anyonlike. As such, the type of interacting topologi- 
cally ordered state with fractional statistic quasiparticles 
such as the fractional quantum Hall (FQH) states, has 
not been found or proposed in an interacting topological 
insulator. The discovery of topologically ordered states 
of matter in the absence of the usual external magnetic 
field applied to the system would be of tremendous in- 
terest to the condensed matter community and could po- 
tentially revolutionize fields such as topological quantum 
computation. 

Interestingly, it has recently been suggested [7HTU] that 
a Chern insulator, i.e., a zero magnetic equivalent of the 
quantum Hall effect [TT], with strong interactions (com- 
pared to the bandwidth) has as a ground state, a frac- 
tional Hall effect at filling 1/3 as it ground state. By 
diagonalizing the Hubbard interaction in a model with 
an almost-flat band dispersion, the authors of [3-^ have 



found, for small sizes of the system, an almost 3-fold 
degenerate ground state, reminiscent of the fractional 
quantum Hall state on a torus. However, the origins 
of such a degenerate state can be either a charge-density 
wave or a fractional quantum Hall state. Arguments were 
given for why the state should be a FQH state, but they 
do not differentiate qualitatively between the FQH state 
and a charge density wave, as will be shown later. The 
very small sizes of the systems considered in the previous 
works and the lack of qualitatively unequivocal calcula- 
tions that prove the existence of a fractional Chern insu- 
lator state make such a state still elusive. In the present 
paper, we present the proof of principle that a z/ = 1/3 
FQH state exists in the nontrivial Chern insulator sub- 
ject to a Hubbard interaction. 

By working in the flat-band limit where the kinetic en- 
ergy is zero and the band gap can be set to infinity, we 
present multiple results that strongly suggest the pres- 
ence of a = 1/3 FQH state. We show that a repulsive 
interaction Hamiltonian has a 3 fold-degenerate ground 
state at filling 1/3 separated from the excited states by 
a gap which we show approaches a finite value in the 
thermodynamic limit. Upon the adiabatic insertion of 
a magnetic flux quantum, the degenerate ground states 
flow into each other - and cross without level repulsion at 
large particle number. The flux insertion and the 3-fold 
ground state degeneracy have already been found in two 
earlier papers, but, without proving the existence of a 
finite gap in the thermodynamic limit, they could very 
well be the hallmark of a charge-density wave (CDW). 
By going to larger sizes than previously achieved, we 
show that the energy gap between the 3-fold degener- 
ate ground state and the first excited states seems to re- 
main finite in the thermodynamic limit when the single- 
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particle Hamiltonian is a nontrivial Chern insulator. In- 
deed, a commensurate CDW at fractional filling would 
also give rise to a gaped phase. Thus, the gap argu- 
ment might be insufficient to prove the existence of a 
Laughlin-like phase. The ground states occur at lattice 
momenta consistent with an emergent generalized Pauli 
principle forbidding the presence of more than 1 particle 
in 3 consecutive orbitals. We compute the momentum 
space density n{k) of each of the 3 degenerate ground 
states and find it to be a constant approaching 1/3 for 
all k, another clue that the ground state is not a CDW 
but a true FQH state. We then show that the quasi- 
hole excitations of these states resemble Laughlin-FQH 
excitations: Their total counting matches exactly that 
of the Laughlin quasihole states on the torus. We also 
present a heuristic Pauli principle which, in some cases, 
counts the number of quasihole states for each momen- 
tum sector. The matching of the excitation properties 
is a clear signal that the state observed is a FQH state: 
Charge density wave states would not have an excita- 
tion spectrum resembling that of a FQH state. We then 
compute the entanglement spectrum (ES) for the ground 
state of the repulsive Hamiltonian, find a large entan- 
glement gap, and show that the counting of the entan- 
glement states below the gap matches that of Laughlin 
quasiholes. The entanglement gap remains finite in the 
thermodynamic limit, and closely tracks the gap in the 
energy spectrum as the system is tuned through a phase 
transition of the single-particle Hamiltonian between the 
nontrivial Chern insulators and a trivial band insulator. 
We also discuss the dependence of the state on the sym- 
metries and aspect ratio of the lattice and on the param- 
eters of the Chern insulators. Both the energy spectrum 
and the entanglement spectrum of the interacting spec- 
trum change fundamentally as the single-particle Chern 
insulator Hamiltonian undergoes a phase transition and 
becomes a trivial band insulator. In the end, we discuss 
the analytical principle behind the counting of quasihole 
states on the lattice. 



THE MODEL AND ITS SYMMETRIES 

Several Chern insulator models with an approximate 
flat band have been recently proposed [SI HOI [12] ■ For our 
model, we pick the Chern insulator on a checkerboard 
lattice, first introduced in [3 [12] • This model already 
exhibits weak dispersion of the bands, but because we 
work in the flat-band limit, this is not essential to our 
calculation - we could have picked an alternate single- 
particle Hamiltonian. As written in [T^j the one-body 
Hamiltonian reads Hi = ^^^(cIa' 4B)'^i(^)(cfeA, Cfcs)^ 
with A, B being the two sites in the unit cell. 
The Hamiltonian matrix can be expressed in terms 
of the 3 Pauli matrices /ii(fc) = '}2ndi{k)ai where 
the di{kys are dx{k) = 4ii cos(0) cos(fca;/2) cos(fcy/2). 
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FIG. 1. Low energy spectrum for A'^ = 6, 10 and 12, A'^^ = 
N/2,Ny = 6. The energies are shifted by Ei, the lowest 
energy for each system size. We only show the lowest energy 
per momentum sectors in addition to the 3-fold ground state. 
We note the good ground state degeneracy, even for relatively 
small system sizes such as 6 particles 



dy{k) — 4<i sin(0) sin(fca;/2) sin(fcy/2), d^ = 2t2{cos{kx) — 
cos(fcj,)) -I- M In the original model [12] there is an ad- 
ditional diagonal term -At^ cos{kx)cos{ky) - which is of 
course relevant for the energy (and shrinks the disper- 
sion of the bands, thereby making them flatter) but does 
not matter for the eigenstates. Since we are diagonal- 
izing in the flat-band limit, we neglect this term. 
is the phase factor added to the nearest neighbor hop- 
pings, while the parameter M is a mass added in order 
to drive the transition from a topological Chern insula- 
tor (for M = 0) to a trivial atomic limit insulator when 
M — > ±oo. The model always has an energy gap (for 
ti,t2,(j) not vanishing) with the exception of the points 
kx = 0, ky = TT, M = -4*2 and kx = ky = 0, M = 4*2 
which are gapless and where the phase transitions be- 
tween the atomic limits M — >■ ±oo and the Chern insula- 
tor phase occur. For \M\ < 4|i2| the filled valence band 
has a Hall conductance of 1. The single particle Hamilto- 
nian matrix has the following symmetries: inversion with 
identity inversion matrix hi{—kx,—ky) = hi{kx 



ky), 



well as (at M = 0) a certain type-of particle hole symme- 
try coupled with a C4 rotation and a mirror operation: 
<Tzhi{kx,ky)az = —hi{ky,kx)- Unfortunately, due to the 
presence effractions k/2, the model in [T2] is not in Bloch 
form. To render it in Bloch form, we perform the gauge 
transformation Cks — ^ CkB exp(— ^(fcj, — ky)/2) to obtain: 



h2{k) = 



^ii(fc) hi2{k) 
hUk) -hn{k) 

hi2{k) = tie"f'{l + e^C'y^''-^) + tie~'"f'{e'''y + e 

hii{k) = 2t2{cos{kx) — cos(fcj,)) + M 



(1) 



The inversion symmetry of hi{k) translates into an- 
other symmetry of h2{k) given by {k)h2{k)U{k) ~ 
h2{—k) with U{k) a diagonal 2x2 unitary matrix with 
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g-«(fex-fcH)/2 on the diagonal. We now fractionally 
fill the valence band of this insulator and add interac- 
tions. The existence of a FQH state in a Chern insu- 
lator should not be taken for granted. Unlike the Lan- 
dau level, a band insulator has nonzero bandwidth and 
a Berry phase distribution which cannot be made uni- 
form over the full Brillouin zone. Fractionally filling the 
band will, even in the presence of large interactions, al- 
low the particles to cover only part of the Brillouin zone 
and hence feel only part of the Berry curvature. To 
eliminate the effect of the band curvature, and to allow 
the filled particles to democratically sample the whole 
Brillouin zone, we always work in the flat-band limit 
of a topological insulator. This corresponds to keep- 
ing the single particle eigenstates of h2{k) but putting 
the energies of the occupied bands to be an arbitrary 
energy zLEq where Eq > 0. At the Hamiltonian level, 
we transform from /i2(fc) = E^{k)P^{k) + E+{k)P+{k) 
to h^^{k) = -EoP-{k) + EoP+{k) where P± are the 
projectors onto the occupied and unoccupied bands. As 
such, the energy difference between the valence and con- 
duction bands can be made large without changing the 
eigenstates of the system. We diagonalize the interaction 
Hamiltonian directly in the filled band, neglecting the 
conduction band. This is similar to the Lowest-Landau 
Level (LLL) projection in the usual fractional auantum 
Hall effect. A nice feature of the checkerboard lattice 
model is that the Hubbard interaction is fixed and has a 
simple form: 



-f^intcraction — ^ ^ Tl-iTLj 
<i]> 



(2) 



where «, j are nearest neighbor sites. One can see that 
the interaction only couples A with B sites. This is not 
the case for multi-orbital models, in which case there is 
both an on-site and nearest neighbour interaction. Upon 
Fourrier transform, and gauge transformation, the inter- 
action reads: 



- V J- ~ 



N 



where 



Vk.MMM = (1 + e'('=^''-'=^''))(l + e-^C^--'^^^)) (3) 

In passing, we remark that the present interaction and 
matrix elements are, in principle, very different from the 
pseudopotential Hamiltonians that give rise to the frac- 
tional quantum Hall effect in the lowest Landau level. 



NUMERICAL PROCEDURE AND RESULTS 

We now diagonalize this Hamiltonian for N particles 
in a Nx x Ny lattice, where, for the ground state we have 
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FIG. 2. Evolution of the 3-fold degenerate ground state upon 
flux insertion along the y direction at A*' = 10, A'^^ = 5 and 
Ny — 6. The 3-fold degenerate ground states spectral flow 
into each other (inset) separated at each point in the flux in- 
sertion from the first excited state (the only one of the excited 
states shown here) which does not exhibit spectral flow with 
any of the other states. 



Nx ■ Ny = 37V (we concentrate on the = 1/3 filling), 
while for the quasihole excitations we have N^-Ny > 3N. 
This setup is quite different from the fractional quan- 
tum Hall effect on a lattice P"5HT5] where the lattice size 
can be changed while keeping fixed both N and v. In 
our case, only the lattice aspect ratio might be tuned 
in some cases (not all of them are accessible for a given 
N and v). All the numerical calculations are performed 
with i2 = (2 - \/2)/2ii as discussed in [H]. The total 
translation operators in the x, y directions commute with 
both the single and many-body Hamiltonians and hence 
the eigenstates are indexed by total momentum quan- 
tum numbers {Kj., Ky) which are the sum of the momen- 
tum quantum numbers of each of the N particles modulo 
{Nx.Ny). The basis states are Ojli 7^ ■ • ■ ^1 ^ '^^^ 



(we work in the "LLL", and the 7^ -'s are the cre- 
ation operators for a particle of momentum k in the 
valence band). When acting on the basis states, the 
%c. = ^-,a.n-,fe ^here is the a = A,B com- 

ponent of the eigenstate of the occupied band of h2{k) or 
h2^{k) (they have identical eigenstates). Diagonalizing 
directly in the valence band provides for large numerical 
efficiency. The inversion symmetry of the single particle 
problem is maintained at the level where the many-body 
interaction is taken into account. Thus the spectrum has 
an exact {Kx,Ky) — >■ {—Kx,—Ky) symmetry which can 
be used as checkup. We first perform all the calculations 
for the M = nontrivial Chern insulator, and then drive 
the fractional Chern insulator to a phase transition by 
increasing M . 
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FIG. 3. Energy gap A for different system sizes and aspect 
ratio. The gap is defined as the difference between the energy 
of the first excited state and the highest energy of the 3-fold 
ground state manifold. In each case, A^^; = 3N/Ny. 



FIG. 4. Ratio between the energy spread 5 of the 3-fold 
ground state manifold and the energy gap A for various sys- 
tem sizes and aspect ratio. The spread is defined as the dif- 
ference between the highest energy and the lowest energy of 
the 3-fold ground state manifold. In each case, A^^, — 3N/Ny. 



Degenerate Ground State and its Properties 

In Fig. [T] we show the spectrum of the system for sev- 
eral sizes N = 6,10,12 and aspect ratios {Nx,Ny) — 
(3, 6), (5, 6), (6, 6) in the M = nontrivial Chern insula- 
tor case . The choice of aspect ratios will be explained 
later. We observe a 3-fold degenerate ground state at 
lattice momenta {(0, 3), (0, 3), (0, 3)}, {(0, 1)(0, 3)(0, 5)}, 
{(0,0), (0,0)(0,0)} respectively. The principle that de- 
termines these momenta will be explained later. The 
insertion of flux quantum in either of the x, y directions 
can be performed by letting each single particle momen- 
tum k^ y — >■ k^ y + ^/N^y wlth ^ S [0,27r]. For the case 
N — 10, Fig.|2]shows that the 3 degenerate states (which 
occur at the same lattice momentum and are hence split 
by the interaction) experience spectral flow into each 
other (as also previously in Ref. [8]). In the presence 
of a gap, 3-fold degenerate ground states that experience 
spectral flow in a = 1/3 fractionally filled Chern insula- 
tor are enough to guarantee the presence of a Uxy = 1/3 
fractional quantum Hall state. Upon the insertion of 3 
fluxes, the ground states move back to the initial config- 
uration, but, as the filled band has Chern number unity, 
one electron has been transferred across the gap, giving 
a Hall conductance of axy = 1/3. These arguments have 
already been presented in [3 |B] but they are valid only 
in the presence of a spectral gap, which was not proved 
in [3 H] • In the absence of such a proof, the effects pre- 
sented in [71 IS] could also occur in a 1/3 charge density 
wave. The existence of a gap in the thermodynamic limit 
is complicated by commensuration effects on the lattice. 
Because of the system sizes that can be accessed numer- 
ically, we cannot give a quantitative value for the size 
of the gap in the thermodynamic limit. Our results do, 
however, strongly suggest that this value should be finite. 

The finite-size scaling of the gaps A is presented in 



Fig.[3j Several crucial trends are visible. For Ny = 3, the 
gap is finite but decreasing as we increase N (or equiva- 
lently — 3N/Ny). Its decrease does not mean that the 
FQH state is compressible in the thermodynamic limit. 
Indeed, the thermodynamic limit here does not corre- 
spond to a 2-dimensional system: for — > oo, if Ny is 
kept fixed at 3, we reach the case of a one-dimensional 
system. This system, even at single particle level, cannot 
have a Chern number- as discussed later, Chern numbers 
appear only in the true 2-dimensional limit. As such, 
it is expected that for fixed A^j^, the gap diminishes in 
the thermodynamic limit. Going to Ny = 6 shows a 
large increase in the gap from Ny = 3 exactly because 
the aspect ratio changed and the system is more two- 
dimensional. However, if we keep Ny fixed, the gap will 
also start diminishing as we go to the thermodynamic 
limit, a fact clearly shown in Fig. [3j When we go to 
Ny — 9, the gap increases again. As such, we are con- 
fident that the gap remains open and scales to a finite 
value for N^/Ny — >■ finite, N ^ oo 

The energy spread is plotted in Fig. |4j If £'i , i?2 , E3 
are the energies of the 3 quasi-degenerate ground states 
and E4 is the energy of the first state above the lowest 3 
states, then the spread S — ~ Ei while the gap A = 
E4 — E^. Fig. |4] reveals the degree of degeneracy relative 
to the size of the gap as a function of A^ = N^ ■ Ny/3 for 
fixed Ny . We see that the degeneracy gets better and 
better for larger sizes. 

We have also plotted the momentum-space density of 
the 3 degenerate states in the Fractional Chern insulator. 
In a FQH state on the sphere or on the torus, the occupa- 
tion of every angular momentum orbital is constant and 
equal to 1/3 (the filling factor) in the thermodynamic 
limit. In the current case, we expect in finite size devia- 
tions from the constant density 1/3 scenario, due to the 
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FIG. 5. Occupation number of each of the single-particle 
momentum orbitals for the three A'^ = 10 ground states (which 
occur at different total momenta) at A^^; — 5, Ny — 6 . Notice 
that the occupation number is uniform and close to the filling 
factor 1/3 



fact that the Berry curvature is not uniform in the Bril- 
louin zone. We however, also expect that the deviations 
not be very large because, in the end, the state must be 
an incompressible liquid. The standard deviation of the 
momentum space density strongly depends on the aspect 
ratio, typically ~ 0.1 for Ny — 3 and ^0.02 for Ny — 6 
(almost independent from N). The results are presented 
in Fig. [sj We observe that the occupation numbers of 
each of the momenta of the 3 degenerate ground states 
for the TV = 10 particle state are close to 1/3, as ex- 
pected for an incompressible liquid at this filling. As the 
3 ground states occur at different momenta, the fact that 
each of the ground states has a uniform n{k) close to 1/3 
for any value of k removes the possibility that these states 
are charge-density waves. 



Quasihole excitations 

To iron out any doubt about the existence of a FQH 
state, we now investigate the quasihole excitations of the 
system. The excitations of a FQH state are very specific 
in the sense that they carry fractional 1/3 charge and 
exhibit fractional 1/3 statistics. This would not happen 
in a charge-density wave state. The quasihole excita- 
tions arise by inserting a flux in the FQH state. These 
properties endow the excitations with a state counting 
specific to the FQH universality class at a certain filling 
independent of the specific Hamiltonian to be diagonal- 
ized or even of the form of the model wavefunction. In 
the current case, it is expected that the actual form of 
the ground state wavefunction is not of a simple known 
Laughlin form. However, it is also expected that the FQH 
at filling 1/3 is the Abelian FQH 1/3 state, whose uni- 
versal properties (after integrating out the fermions) are 
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FIG. 6. Low energy spectrum for A'^ = 9, Nx = 5, Ny 
The number of states below the dashed line is 19 in sectors 
where Ky mod 3 = 0, 18 in other sectors, as expected from 
the analytical results. The total counting matches the one 
expected for Laughlin quasihole states. The counting for each 
momentum sector is given by the generalized Pauli principle 
and two dimensional to one dimensional unfolding presented 
in the next section. 



described by the topological Abelian Chern-Simons field 
theory. As such, its universal quasihole counting (even in 
finite size) should be identical to that of the usual Laugh- 
lin quasihole states. In fact, as our interaction represents 
a generic Hamiltonian (similar to the Coulomb Hamilto- 
nian in the FQH case), all we can hope for is that the 
set of quasiholes is separated by a finite gap from higher 
energy non-universal excitations that are not described 
by the topological field theory. 

The FQH i/ = 1/3 quasihole counting is identical to 
that of particles with Haldane 1/3 statistics in both the 
thermodynamic limit and in the finite size systems. This 
counting can be directly obtained from the topological 
Chern-Simons field theory and the assumption of bulk- 
edge correspondence: In a system with the boundary, 
the Abelian Chern-Simons field theory is not gauge- 
invariant - it is missing an edge piece which is a boson 
at compactification radius \/\Jv. The modes of this bo- 
son are the edge modes, a gapless conformal field theory 
whose excitations, when placed in a finite-size box exhibit 
Haldane Xjv statistics: Every particle must be separated 
from another particle by at least Xjv orbitals. This is 
called (1, Xjv) generalized Pauli principle, which will be 
presented in the next section. However, the bulk-edge 
correspondence renders this counting also the same as 
that of the quasiholes in the bulk of the FQH state, and 
represents a hallmark of the = 1/3 FQH state. 

We now compare such a counting with that obtained 
from numerics. For the model FQH Hamiltonians such as 
the Haldane pseudopotentials, the quasiholes are just the 
degenerate ground states of the system when fluxes are 
added to the system. In this sense, the ground state is 
just the highest density state of a set of other degenerate 
states which occur at higher fluxes and which describe 
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the quasiholes. As we cannot add unit fluxes to a ground 
state with a set aspect ratio and still keep the transla- 
tional symmetry intact, we choose to analyze quasiholes 



by either increasing — >■ + 1 (or Ny 



1) or 



by removing particles from a ground state and keeping 
the aspect ratio fixed - hence our system sizes contain 
at least 3 fluxes above the ground state. A generic ex- 
ample is that shown in Fig. [6] in which = 9 particles 
reside in x Ny = 5x6 orbitals, 3 fluxes more than the 
ground state at filling 1/3. We diagonalize the interact- 
ing Hamiltonian for this number of particles and observe 
(Fig. [6]) that the spectrum splits into two parts sepa- 
rated by a clearly visible and unambiguous energy gap. 
The counting of states below the gap equals the counting 
of quasiholes of a 1/3 fractional quantum Hall state of 
9 particles in 30 orbitals (explained in the next section). 
Even more, the counting of states per momentum sector 
also matches that of the generalized Pauli principle for 
the quasihole states ^ previously obtained in the FQH 
context. This will also be explained in the next section. 

While for N = 9, Nx X Ny = 5 X 6 the gap that sep- 
arates the quasihole manifold from the spurious states 
is clear, this is not always the case. In Fig. [7| we plot 
the spectrum of the N = 7, N^ x Ny = 8 x 3 problem: 
for most momentum sectors we find a clear gap, visible 
by eye, between a low-energy quasi-hole manifold and 
that of spurious, high-energy states. However, for the 
momenta (K^jKy) = (2,0) (and for its inversion sym- 
metric (6, 0)) we see one state (at energy around 0.21) 
which lies in the gap. It is a-priori unclear whether, just 
based on the numerical information, this state belongs 
to the quasihole manifold or not (the analytic counting 
obtained in the next section would indeed tell us that 
it does belong to the quasihole manifold). To find this 
out, we look at the level spectral flow of the states at 
{Kx,Ky) = (2, 0) upon the flux insertion kx kx + ^/N^ 
in Fig. [Sj where we insert flux in the x direction. We 



note that the state at energy 0.21 exhibits level spec- 
tral flow with the states below it, and hence it belongs 
to the quasihole subspace. Note that the first two ex- 
cited states (around energy 0.24) do not mix with the 
quasihole manifold upon flux insertion and even exhibit 
level repulsion between themselves. For example, when 
$ = 7r/8 the lower energy manifold is clearly separated 
from the spurious, high-energy one. The flux insertion 
tells us that the state which, in the periodic boundary 
condition case was at energy 0.21 indeed belongs to the 
quasihole manifold. 

A similar analysis has been made on several other sys- 
tem sizes and aspect ratios for quasiholes close to the 
1 /3 fllling. In the absence of a pseudo-potential Hamil- 
tonian and since we work on a lattice, several factors 
must be taken into account when looking at quasiholes. 
First, as we are diagonalizing a generic Hubbard- type 
Hamiltonian, we must make sure that we are still close 
to filling 1/3. Taking out too many particles (creating 
too many quasiholes) could potentially lead us to another 
filling factor and diminish or even erase the gap between 
the low-energy quasiholes and the spurious high energy 
states. Since we work on the lattice, care should be taken 
that the single-particle system is not one-dimensional - 
the aspect ration should be roughly balanced. This will 
be discussed later. In most cases when these two condi- 
tions are satisfied, the gap between the quasihole man- 
ifold and the upper spurious states is visible. The to- 
tal number of states (in all momentum sectors) below 
the spectral gap matches exactly the number of states of 
the Laughlin quasiholes in the regular FQH effect on the 
torus of N particles and Nx ■ Ny number of orbitals. This 
is the number of (1, 3)-admissible partitions on the torus 
(as defined below, see Figs. |9] and 10 ) that we can write 



down out of N particles in Nx ■ Ny orbitals. For N 
and Nx, Ny = 5, 6 , respectively, we obtain a total of 6435 
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states while for = 9 we obtain 550 states by counting 
these partitions. This number is identical to the numer- 
ical data (see Fig. [6] and [7| . We have checked a series of 
sizes up to = 10 and have confirmed in each case the 
existence of a large gap between a high energy nonuniver- 
sal sector and a low-energy quasihole sector whose count- 
ing matches exactly that of the Laughlin quasiholes ob- 
tained by counting (1, 3)-admissible partitions. All these 
quasihole states exhibit spectral flow within themselves. 
This proves that the observed state is an Abelian frac- 
tional quantum Hall state with a^y = 1/3, as it has ex- 
citations obeying Haldane 1/3 statistics and having the 
counting of a U{1) boson. The overall quasihole count- 
ing matches the (1,3) generalized Pauli principle. But in 
most cases, this statement is also valid for the counting 
per each momentum sector. The (1,3) generalized Pauli 
principle is generic of the Laughlin FQH state at filling 
1/3 and will be presented below. 



HEURISTIC COUNTING RULE: PAULI 
PRINCIPLE AND CRYSTAL MOMENTA 

We now present a counting rule for the total lattice 
momentum at which the degenerate ground states occur. 
A more detailed version will be developed in an upcoming 
paper [T7] . The 2-dimensional lattice of sites x Ny with 
kx € [0, Nx — 1] and ky e [0, Ny — 1] can be folded into a 1- 
dimensional lattice of orbitals of momentum A = N^ky + 
kx by placing the x momenta range one after each other 
for the different values of the y-momenta, as per Fig. |9] 
This then represents a map between the momentum of 
the i'th particle orbital kx^i,ky_i and the 1-dimensional 
momentum of the i'th orbital Xi E [0, ■ Ny — 1] . The 
Ai's then form a partition A = [Ai, . . . Aat] which we order 
from highest to lowest. 

This partition is then similar to the orbital momenta of 
the particles of the usual FQH on the torus. On such ge- 
ometry and in the lowest Landau level, the single-particle 
orbitals can be written as 



,FQH 



{\+mN^){x+iy) 



(4) 



xe be 



In the above, we have picked the Landau gauge A = 
BxCy. lb is the magnetic length, N^ is the number of flux 
quanta and Ly the system size along y. With this gauge 
choice, A is the orbital momentum along the y direction 
and is such that A € [0, A^^ - 1]. Setting N^ = N^ ■ Ny, 
we can formally relate the FQH momentum A and the 
linearized momentum of the FCI A = N^ky + k^. 

Using this analogy between FQH on the torus and FCI, 
we now index each many-body state by a single partition 
(called "root" partition [IB]) that satisfies a generalized 
Pauli principle (also called an "admissibility condition") 
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FIG. 9. Momentum of the 3-fold degenerate ground state 
for the fractional Chern insulator at filling u — 1/3. The 2- 
dimensional lattice of momenta kx,ky with kx G [0,A'^i; — 1] 
and ky G [0, A^^ — 1] (in units of 2n/Nx,2n /Ny respectively) 
where we have taken to be a multiple of 3 (since ■ Ny = 
3N, this can always be done) is unfolded in a 1-d lattice. The 
total momentum of the ground states is then the same as the 
momentum of the three (1, 3)-admissible partitions possible 
for the number of orbitals N^ ■ Ny and number of particles 
N = Nx- Ny/3 



which does not allow the existence of more than 1 particle 
in 3 consecutive orbitals. In partition notation, this reads 
Xi — Ai+i > 3. Because of the periodicity of the torus, we 
must also make sure that the first and last particles are 
separated by at least 3 orbitals, which reads: Xn + N^ ■ 
Ny-l-Xi > 3. At the ground state filling A^ = N^-Ny/S 
(we pick N^ a multiple of 3) there are only three such 
partitions that correspond to the occupation numbers 
100100100... 100100 {Kx,Ky ^ N^iN^ - 3)/6, N{Ny - 
l)/2 mod {Nx,Ny)), 010010010... 010010 {Kx,Ky = 
Nx{Nx - 3)/6 + N,N{Ny - l)/2 mod (Nx^Ny)), and 
00100100100... 001001 {K^,Ky = N^{N^ - 3)/6 + 
2N,N{Ny - l)/2 mod {N^,Ny)), of the orbitals of mo- 
mentum Ai, A2, . . . Aat. As an empirical observation, it 
was noted in T] that the degenerate ground states ap- 
pear at momenta related by a translation with the same 
number in both x,y directions. This is an immediate 
corollary of the counting principle presented above. In 
the present model, due to inversion symmetry, at least 
two of the ground states are exactly degenerate. The 
principle for finding the ground state momentum reminds 
one of the root partitions used in the usual FQH effect 
[TB] or in the thin-torus limit of the FQH states [T5| . 

At this point, a discussion of the dependence of the 
physics on the aspect ratio of the problem is necessary. 
In cases of finite-size systems, the exact diagonaliza- 
tion data for the momentum of the 3 quasi-degenerate 
ground states matches the momentum counting pre- 
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sented above for aspect ratios for which the sample is 
two-dimensional. The Chern insulator problem has the 
particular property that the filled band only has a Chern 
number (Hall conductance) equal to unity in the case of 
Nx/Ny — > finite. For example, if we let Ny = 1 and 
Nx — > oo, the Chern number of this Hamiltonian will 
not be equal to unity. The Chern number is equal to 
/ dkxdkyTr (P[9fc^P, dk^P]) /(27r)^ and hence Ny cannot 
take small values such as 1, 2 in order for the derivatives 
of the Chern number to be well defined. The Brillouin 
zone mesh of momenta must truly be two-dimensional in 
order for the filled band to have an integer Chern num- 
ber. We hence expect that the interacting problem also 
be sensitive to the aspect ratio when one of the dimen- 
sions becomes much smaller than the other (of course in 
the thermodynamic limit, as long as Nx/Ny remains fi- 
nite, we expect the degeneracy of the FQH states to be 
independent of the aspect ratio). In fact, the interacting 
problem has also a different dependence on the aspect 
ratio. To see this, one can look at the extreme case of 
A^j^ = 1, Nx — 3N. In this case, we are indeed solving a 
one-dimensional problem of orbitals kx = 0, . . . , 3A^ — 1 - 
extremely similar to that of the LLL orbitals on the torus 
- but with interacting Hamiltonian that comes from a 2D 
Hubbard interaction. This Hamiltonian is much different 
from the Haldane pseudopotential Hamiltonian which in- 
volves Clebsch- Gordon coefficients on the sphere and Ja- 
cobi theta functions on the torus. As such, the problem is 
physically distinct at the interacting level, and for these 
skewed aspect ratios we do not expect FQH states in the 
Chern insulator problem. Indeed, in the finite size di- 
agonalizations performed below, we will always aim to 
choose the most symmetric aspect ratio available. As 
observed in Fig. |3] and |4j the cleanness of the results will 
depend on this choice, as it is natural for these finite 
sizes. 

The counting of quasiholes per momentum sector is 
slightly more complicated, and we have not found the 
full counting rule per momentum sector. The map be- 
tween the two-dimensional momentum numbers and the 
one-dimensional orbital partition that allowed us to find 
the momentum of the ground states cannot always work 
for the quasihole states. It sometimes violates the in- 
version symmetry (see the Appendix for such cases) and 
thus it sometimes does not match the numerical results. 
Nevertheless, the total number of quasiholes, irrespective 
of the momentum sectors, always matches the counting 
of partitions. This matching of the exact diagonalization 
states and the 1.3 - admissible partitions is a rather re- 
markable demonstration of the = 1/3 character of the 
states. Even more interestingly, the Pauli principle and 
the map between the two-dimensional momentum num- 
bers and the one-dimensional orbital partition works for 
each momentum sector (not only for the total number of 
quasiholes) in the great majority of the cases studied. As 
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FIG. 10. Example of total momentum counting for the 3- 
fold degenerate ground states of the N = 4 (top) and = 6 
(bottom) problem. Every time the number of particles A'^ is a 
multiple of both kx,ky, the 3-degenerate ground states occur 
at the same total momentum and are expected to be split 
by the interaction for finite-size samples. Notice that for the 
TV = 4 problem, the total momenta at which the ground state 
occurs are related by inversion symmetry, as (1, 2) = (—2, —2) 
mod (3,4) 
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FIG. 11. Example of total momentum counting for the quasi- 
holes of the = 2 7V^ = iVy = 3 problem. 



counting principle at work for N — 2, Nx 



N, 



an example, we present in Fig. 11 the generalized Pauli 



ENTANGLEMENT SPECTRUM 

We now turn back to the 3-fold degenerate ground 
state of the system. There are several attempts to ob- 
tain an analytical expression of the Laughlin state for 
the fractional Chern insulators [ini HHj ■ Unfortunately, 
none of them can be used in our case, so we cannot per- 
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form any wavefunction overlap calculations to compare 
their analytical expression with our system ground state. 
Nevertheless, we can still show that the ground state con- 
tains by itself, information about the Abelian fractional 
1 /3 character of the excitation spectrum. This is actually 
a far better probe than an overlap, since the 1/3 charac- 
ter is the feature, not the analytical expression. To do 
this, we use the recently developed entanglement spec- 
trum [m [22] which for a single nondegenerated ground 
state |\[') can be defined through the Schmidt decomposi- 
tion of l^*) in two regions A, B (not necessarily spatial): 



|*> = E' 



(5) 



where |*/) = (*f l^-f ) = The exp(-eO and 
) are the eigenvalues and eigenstates of the reduced 
density matrix, pA — T^t^bP, where p — \'^) (^'| is the 
total density matrix. There is no generalization of ([5| 
to degenerate ground states. Still the definition of the 
entanglement spectrum through the reduced density ma- 
trix can be extended to that case. While several schemes 
can be proposed, it has been observed in [55] that the 
incoherent summation over the degenerate ground states 
P — X)i l^i) is good candidate for this generaliza- 
tion. This combination builds a density matrix which 
commutes with the total translation operators which is a 
desired feature to sort the with respect to the momen- 
tum quantum numbers. 

Depending on the space where the system is split into 
A and B, be it real, momentum or particle space, dif- 
ferent aspects of the system excitations will be revealed 
through the ES. It was proven that if the regions A, B are 
regions of particles [35], the particle entanglement spec- 
trum hence obtained by tracing over the positions of a 
set of B particles gives information about the number of 
quasiholes of the system of Na particles and the number 
of orbitals identical to that of the untraced system. In the 
case of the usual FQH, the particle entanglement spec- 
trum of a model state contains an identical number of 
levels (i.e., the number of non zero eigenvalues) as those 
of the quasiholes. Thus the counting of non zero eigen- 
values does not suffer from the sometimes uncontrolled 
finite-size effects, as observed with state decomposition 
in the momentum space |211 123| . For this reason we will 
use state decomposition in the particle space as a probe 
for the fractional Chern insulator. 

Away from the model states, like the Coulomb ground 
state, the ES may exhibits an entanglement gap [2T1 [21] . 
It separates a low-energy structure with perfect quasihole 
counting and and a high entanglement energy nonuniver- 
sal part. But a clear and significant gap is not always 
observed, even for the = 1/3 Coulomb state. 

For the fractional Chern insulator, the situation is sur- 
prisingly much better: We observe a clear, large entan- 
glement gap between low entanglement energy levels and 
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FIG. 12. Particle entanglement spectrum for A'^ = 10, Nx 



5,iV„ 



6 and Na 



5. The number of states below the 



dashed line is 776 in all the kx = 0, and 775 in the other 
sectors. This is in agreement with the (1,3) Pauli principle 
for every momentum sector. 



the high entanglement energy levels like those observed 
Moreover, the counting of the levels 



in Figs n2^ and 13 



below the gap is identical to the counting of quasiholes of 
Na particles in A^^, • Ny orbitals. In both examples shown 
in Figs. [T2| and |13[ it matches per each momentum sec- 
tor that of (1, 3)-admissible partitions of A'^ particles in 
Nx ■ Ny orbitals. We find it very revealing that the frac- 
tional Chern insulator ground states obtained here con- 
tain much clearer information (large, clear entanglement 
gap) than the ground states of the Coulomb interaction 
in the FQH. The entanglement spectrum shows that the 
ground states by themselves contain essential information 
on the fractional nature of the excitations in the frac- 
tional Chern insulator. The current clean application of 
the entanglement spectrum also shows that this quantity 
is fundamentally useful toward revealing the physics of 
strongly-correlated states besides the usual FQH model 
wavefunctions for FQH states. 

The ES also provides some insight about the system 
when the number of sites in one direction is equal to 3. 
A clear entanglement gap is observed and the counting 
below this gap matches the expected one but only for 
Na = 2 and A^^i = 3. For larger values of A^^, the num- 
ber of levels below the gap is lower than expected. This 
strongly suggest that for the small aspect ratios where 
the ration Ny/N^ — 0, the state is not a fully developed 
fractional Chern insulator - or indeed, by our previous 
arguments, the single particle problem is not even a well- 
developed integer Chern insulator. In the case of the 
other aspect ratios, all the ESs for any value of Na and 
any system sizes up to A^ = 12 perfectly match the pre- 
dicted counting. For the cases such as N^ x Ny = 6x5 and 
6x6, the (1, 3) Pauli principle counting matches the ES 
data for each momentum sector for all A^^ — 1, 2, 3, 4, 5. 
The perfect match clearly shows that this state is a Frac- 
tional Chern Insulator. 
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FIG. 13. Particle entanglement spectrum for A'^ = 12, A^^ = 
Ny = 6 and Na = 4. The ^'s are the entanglement energies. 
The number of states below the dashed line is 741 in mo- 
mentum sectors where mod 2 = A^j, mod 2 = and 728 
elsewhere. The total number below this line (26325) exactly 
matches the one predicted by the counting rule. 



TRANSITION TO THE TRIVIAL INSULATOR 



Tuning the mass term above the 4i2 threshold yields 
an insulator topologically equivalent to the atomic limit, 
with zero Hall conductance. Partially filling this insula- 
tor at say M = 6t2, we find no clear sign of a 3-fold de- 
generate ground state, and no clear gap to the excitation 
spectrum (see Fig. 14). In the interacting problem, we 



would think that the gap re-arrangement happens more 
quickly than M = At2- This hypothesis is natural since 
the many-body gap is expected to be more sensitive than 
the single-particle one. The evolution of the many-body 
gap with the parameter M is plotted in Fig. [15] In the 
atomic limit M — >■ cc, B sites have an energy —M and all 
the particles are strongly localized on those sites. If all 
Bsites were occupied and A sites unoccupied, the filling 
would be 1/2, and hence as we are at a smaller filling 1/3, 
and as the interaction couples sites A and -B, the ground 
state will be highly degenerate. We see that the many- 
body gap as a function of M has a collapse at exactly the 
value of M at which the one-body transition takes place. 
We also plotted the entanglement gap as a function of M 
in Fig. [TBj The entanglement gap is defined as the min- 
imum separation over the different momentum sectors, 
between higher entanglement energy level belonging to 
the Laughlin counting and the lower entanglement en- 
ergy level that does not belong to this counting. PES 
tracks the energy gap and its collapse also at the same 
value as the many-body and single-particle energy gaps. 



FIG. 14. Low energy spectrum for different values of M at 
A'' = 10, Nx ~ 5 and A^j, = 6. We show only the two lowest 
energy per momentum sector. All energies are shifted such 
that the ground state energy is for each M value. 
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FIG. 15. Energy gap A as a function of M at filling = 1/3 

for Af = 8 and A^ = 10 with A^^ = N/2, Ny = 6. For both 
system sizes, the transition occurs at M — 4t2. 
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ANALYTICAL COUNTING OF QUASIHOLES 
PER MOMENTUM SECTOR 



The heuristic generalized Pauh principle counting of 
the previous section gives the correct (matching the 
data) total counting of quasiholes states manifold at 
any N, , Ny (as long as the aspect ratio is two- 
dimensional) . In most cases, although not in all, it also 
gives the correct counting per each momentum sector. 
The generalized Pauli principle is based on the corre- 
spondence between a 2D momenta map to ID orbitals, 
which is heuristic. In this section we provide some ana- 
lytic results of the counting per momentum sector that 
has been tested to work for all the momentum sectors 
and is not based on a 2-d to 1-d folding. The counting of 
the total number of quasihole states of a 1/3 state of N 
particles in ■ Ny orbitals on the torus can be obtained 
by counting the number of ways of putting N particles in 
Nx ■ Ny orbitals with the restriction that there should be 
no more than 2 particles in 3 consecutive orbitals ( (1,3) 
admissible partitions-generalized Pauli principle) on the 
sphere: 



#qh sphc 



N- 



Nj:Ny + 2-'3N (6) 



where n is the number of quasiholes, and then subtract- 
ing the configurations that violate the (1,3) generalized 
Pauli principle once we make the system periodic to go to 
the torus. These configurations come in only three pos- 
sible ways and can be written in terms of the occupation 
number of orbitals 

0100(all (l,3)-admissible configurations)001 
100(all (l,3)-admissible configurations)0010 
100(all (l,3)-admissible configurations)001 (7) 

These represent the only sets configurations that are al- 
lowed on the sphere but that have more than one particle 
in 3 consecutive orbitals if the orbital space is made peri- 
odic. The number of the first two are equal by inversion 
symmetry and reads: 



N 



-2 
n' 



while the third one reads 

' N ~2 + n'^ ■ 
n" 

The result then reads: 



n' = N^Ny + l-2,N (8) 



N^Ny + 2-m (9) 



_ {N,Ny-2N-l)\ 
#qh torus - i^.i^v m{N,Ny - 37V)! 



(10) 



This is the total number of quasiholes of a 1/3 state 
on the torus In a translationally invariant system, and 



for a FQH incompressible and featureless, topological 
state, there is no particular reason why some momenta 
should have an occupation number different from oth- 
ers. Indeed, both the interacting-particle and the single- 
particle models are translationally invariant, and as such 
in the featureless liquid ground state, momenta occupa- 
tion numbers should be the same. Hence, on a finite-size 
lattice for which we have N^ ■ Ny momenta, we would 



expect to have 



(N^Ny-2N-iy. 



states for each momenta. 



NliN^Ny-SNy. 

While this is always certainly true in the thermodynamic 
limit (as if the heuristic Pauli principle), there are cases 



in finite size when 



(Af^Ar„-2Af-l)! 



jV!(jv.Ar„-3jV)! not an integer. This 
means that in finite size cases, commensuration effects 
weakly modify the finite size counting of each sector so 
that the requirement that there be an integer number of 
quasihole states for each momentum sector is satisfied. 
This is the reason for the finite size changes of the num- 
ber of states in each momentum sector. 

Physical reasoning leads us further. When N is rela- 
tively prime with both N^ and Ny there can be no com- 
mensuration effects, and the number of quasihole states 
per each momentum sector has to be identical. Indeed, 
in this case, ^^fiN" N^-sNy. integer, is equal to the 

number of states per momentum sector provided by our 
heuristic Pauli principle, and matches the numerical data 
for the number of states per momentum sector. This is 
also the asymptotic value of the number of states per 
momentum sector in the thermodynamic limit. 

Finite size effects become important in the case when 
the greatest common divisors (GCD) of N, N^ or N, Ny 
is larger than 1. Three cases can occur. First, 
GCD{N,Ny) > l,GCD{N,N^) = 1: in this case, the 
number of quasihole states is the same for all the Ky 
total momentum that are divisible by GCD{N, Ny) and 
is different from the number of states at Ky total mo- 
mentum not divisible by GCD{N, Ny). A clear example 
of this situation is in Fig. [6} N ~ 9,Nx = 5,Ny = 6 
and hence GDG{N,Ny) = 3. The momenta {K,^,Ky) 
for which KyTnodi — have 19 states in the quasihole 
subspace, whereas all momenta for which ii'j,mod3 ^ 
have 18 states in the quasihole subspace. This is also 
equivalent to the counting of (1, 3) partitions obeying the 
generalized Pauli principle. The case GCD{N, N^) = 
l,GGD{N,Nx) > 1 is obviously the x ^ y of the 
case presented above. Second when GCD{N, Ny) = 
GGD{N, Nx) > 1, the number of quasihole states is the 
same for all K^, Ky which are simultaneously both divisi- 
ble by GCD{N,Ny){^ GGD{N,Nx)) and different from 
the number of quasihole states at K^ , Ky when either K^ 
or Ky is not divisible by GCD{N, Ny). An example of 
this latter situation is in the counting of the entanglement 
spectrum states of the iV = 12 {N^ = 6, Ny — 6) par- 
ticles ground state and A^^ = 4 in Fig. [13] The counting 
of states in the entanglement spectrum should be iden- 
tical with that of Na = 4 quasiholes in A^^; • A^j, = 36 
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orbitals. As such, GCD{N,N^) = GCD{N,Ny) = 2. 
Momenta {K^,Ky) = (0,0), (2,0), (4,0), (0,2), (2,2), 
(4,2), (0,4), (2,4), (4,4), which are both divisible buy 
GCD{N, iVj.) have 741 quasiholc states, whereas all other 
momenta, in which either or Ky are not divisible by 
GCD{N,Nx) have 728 quasihole states, for a total of 
26325 states, which matches the total number of states 
of the (1,3) Pauli principle. These rules for the count- 
ing of states are physically motivated by the fact that 
momenta on the lattice should be filled democratically 
except in the case when the number of particles exhibits 
some commensuration (expressed by us as the condition 
GCD{N, Nx^y) > 1) with the lattice dimensions. They 
have been checked for a range of particle number (up 
to N = 12) and found to hold in all cases. The rules 
make sense physically: In the case of commensuration, 
the momenta which are divisible by the commensura- 
tion factor (GCD) should all exhibit the same count- 
ing (once the commensuration condition has been estab- 
lished, the commensurate momenta are treated demo- 
cratically). This counting should be different from the 
counting of the momenta incommensurate with the GCD. 
Whenever either GCD{N,N^) or GCD{N,Ny) equal 
unity the Pauli principle counting of states matches the 
numerically found counting. 

The third case, in which GGD{N,Ny) ^ 
GCD{N,Nx) ^ 1, has not been analyzed and compared 
to the numerical data (as the size of the numerical 
computation becomes too large). However, we can 
offer a conjecture for the quasihole state counting per 
momentum. (Kx,Ky) momenta which are not multiples 
of either GCD{N,Ny) nor GCD(N,N^) have identical 
counting, different from all others. The momenta for 
which GCD{N,N^) divides but GCD{N,Ny) does 
not divide Ky all have the same counting, different 
from any other. This is similar for momenta for which 
GCD{N,Ny) divides Ky but GCD{N,N^) does not 
divide K^- Finally, momenta for which K^ is divisible 
by GCD{N,N^) and Ky is divisible by GCD{N,Ny) 
have same counting, different from the rest. 

CONCLUSION 

In conclusion, we have shown that the ground state 

of the v = 1/3 flat-band Chcrn insulator in the pres- 
ence of repulsive interactions is an incompressible state 
with Hall conductance 1/3, and that quasihole excita- 
tions satisfying fractional statistics. The presence of the 
excitation whose counting is identical to those in the FQH 
u = 1/3 state is a clear and up-to-now missing proof that 
the ground state is indeed a fractional Chcrn insulator 
and not a 1/3 charge-density wave state. We have pre- 
sented a mapping between the lattice momenta and the 
torus orbitals and shown that the counting and total mo- 
menta of the ground states and those of the quasiholes 



can be obtained by employing a generalized Pauli princi- 
ple of not more than 2 particles in 3 consecutive orbitals. 
We have then shown that the entanglement spectrum of 
the ground state also has an entanglement gap and that 
the levels below the gap match in counting those of the 
Laughlin quasihole states. 
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APPENDIX I: QUASIHOLES AND THE 
GENERALIZED PAULI PRINCIPLE 

The counting of quasiholes per momentum sector is 
slightly more complicated, and we have not found the full 
counting rule per momentum sector. Although the total 
number of quasiholes always matches the counting of par- 
titions, the map between the two-dimensional momentum 
numbers and the one-dimensional orbital partition that 
allowed us to find the momentum of the ground states 
cannot always work for the quasihole states. This can 
be easily seen as the partition map between the momen- 
tum of the i'th particle kxi,kyi and the one-dimensional 
momentum (1, 3)-admissible partition A = (Ai, . . . , Ajv) 
with \i = Nx- kyi + kxi does not (for quasiholes) have to 
respect the inversion symmetry that the model studied 
has. To illustrate this, take two particles their mo- 
menta {kxi,kyi), {kxj,kyj) and their orbital momentum 
Ai,Aj satisfying the (1,3) admissibility \\i — Xj\ > 3. 
The inversion property acts at the single-particle level 
and transforms k^^ky —kx,—ky. To make the in- 
verted momentum belong to the positive numbers, we 
must add N^^Ny to the x,y momenta except if the mo- 
menta are 0, in wliic;li case we do nothing. With these 
new momenta, we form the inverted partitions A™" , A"" 
and check whether these violate the admissibility rule, 
i.e. whether |A"'" — A""'[ < 3. We then have to analyze 
seven cases separately out of which three cases turn out 
to be problematic: 

1. kxi ^ 0, ky^ ^ 0, kxj ^ 0, kyj 7^ 0, Af ^' = N^Ny + 

7V,-A„ Aj™ = N^Ny+N^-Xj and lA^^-Af^'l = 
|Aj — Ajl > 3 - inversion symmetry preserved. 

2. fc,, = 0, ky, ^ 0, k^j ^ 0, kyj ^ 0, = NxNy " 

A„ Af ^ = N^Ny + Nx- Xj and lAf'" - Af''| = 
+ K — - and hence inversion symmetry is 
not preserved in the cases, kyj = kyi , k^j = Nx — 
2? A^x 1 or kyj = kyi -\- 1, kxj = 1, 2^ 
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0, A™" = N^Ny 



3. kxi 0, kyi ^ 0, kxj 7^ 0, kyj 

A„ ^ N^- Xj and jA™' 

-^a; — Ai + Aj I - and hence inversion symmetry is not 
preserved in the cases kyi — Ny — 1 , fc^ 



a;™ I 



\N,Ny 



1,2; 



4. 0, kyi ^ 0, fca; 



0, ky 



M 7^ '^i 

A-1 - |A,-A,| > 
3 - inversion symmetry preserved. 



N^Ny- 



Xj and I A™- 



5- — 0, kyi — 0, k^j 0, kyj 0, X™^ — Aj — 0, 
Af " = N,Ny+N,-Xj and lAr^-Aj™! = |7Vj;A^,+ 
-^a; + Aj — Aj I > 3 if Nx > 1 - inversion symmetry 



preserved. 



6. kr 



0,k. 



0, 7^ 0, fcyj — 0, XI 



kxj ~ 

not preserved 



Xj and I A*"" 



< 3 



= A, = 0, 
for Xi — 



2,Nx — 1 then inversion symmetry is 



7. kx 



0,k. 



yi 



XJ — 0, kyj 7^ 0, Aj 



A, - 0, 



a;."" = NxNy- 

3 - inversion symmetry preserved. 



-Xj and |Af"- 



-a;™i 



|7V,7Vj,-A,| > 



Despite the above, in most (but not all) of the cases, 
the counting of states per momentum is still given by 
the (1,3) -partitions of the same total momentum Kx = 
J2iLi ^xi and Ky = J2iLi kyi, even though there is no 
inversion symmetry in the partitions (as proven above; 
we stress that the total counting of quasiholes always 
matches that of (1, 3) partitions, as it should for a Laugh- 
lin state; only the 2D to ID mapping proposed in this 
paper sometimes fails for quasihole states). What hap- 
pens in these cases is that momentum sectors related by 
inversion symmetry both contain the same number of 
partitions which under the inversion operation are not 
(1, 3) admissible. So the effect of the 2D to ID mapping 
on the sectors not respecting the inversion symmetry is 
magically canceled. We find that the matching per mo- 
mentum sector of the exact diagonalization states and 
the 1.3 - admissible partitions to be a rather remarkable 
demonstration of the = 1/3 character of the states. 
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